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[l] Motivated by the physical picture of shape-dependent 
air resistance and, consequently, shape-induced differential 
sedimentation of dust particles, we searched for and found 
evidence of dust particle asphericity affecting the evolution 
and distribution of dust-scattered light depolarization ratio 
(5). Specifically, we examined a large data set of 
Cloud-Aerosol Lidar with Orthogonal Polarization 
(CALIOP) observations of Saharan dust from June to 
August 2007. Observing along a typical transatlantic dust 
track, we find that (1) median § is uniformly distributed 
between 2 and 5 km altitudes as the elevated dust leaves the 
west coast of Africa, thereby indicating unifonnly random 
mixing of particle shapes with height; (2) vertical 
homogeneity of median § breaks down during the westward 
transport: between 2 and 5 km 8 increases with altitude and 
this increase becomes more pronounced with westward 
progress; (3) 5 tends to increase at higher altitude (>4km) 
and decrease at lower altitude (<4km) during the westward 
transport. All these features are captured qualitatively by a 
minimal model (two shapes only), suggesting that 
shape-induced differential settling and consequent sorting 
indeed contribute significantly to the observed temporal 
evolution and vertical stratification of dust properties. By 
implicating particle shape as a likely cause of gravitational 
sorting, these results will affect the estimates of radiative 
transfer through Saharan dust layers. Citation: Yang, W., 
A. Marshak, A. B. Kostinski, and T. Vamai (2013), Shape-induced 
gravitational sorting of Saharan dust during transatlantic voyage: 
Evidence from CALIOP lidar depolarization measurements, 
Geophys. Res. Lett., 40, 3281-3286, doi:10.1002/grl.50603. 

1. Introduction 

[ 2 ] Accurate estimation of net radiative influence of dust is 
difficult because of the complex temporal and spatial 
variations in dust properties. Dust optical properties vary 
with the chemical composition [e.g., Sokolik and Toon, 
1999; Petzold et al., 2009] and size of particles, and can 
change as dust interacts with other substances in the air 
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[e.g., Levin et al., 1996; Trochkine et al., 2003]. In 
addition, the scattering and absorbing properties of dust 
are also affected by particle shape [ Kalashnikova and 
Sokolik, 2002; Dubovik et al., 2006; Nousiainen, 2009]. 
For example, the asphericity of particles can enhance 
forward scattering and reduce back scattering and can 
make single scattering albedo and extinction coefficient 
different from their spherical counterparts. Hence, particle 
asphericity needs to be considered in dust radiative 
forcing calculations. To that end, we recall that particle 
shape affects air resistance (drag) as well and can, there- 
fore, influence the vertical distribution of dust properties 
during transatlantic transport. Hence, the radiative and 
fluid-mechanical aspects of dust transport can be linked via 
particle shape and we now proceed to employ this link. 

[ 3 ] Because of higher drag, at a given volume and mass, 
aspherical particles fall more slowly and hence stay aloft 
longer than their spherical counterparts. Note that this is 
particularly pronounced for lower Reynolds numbers 
(laminar regime). Theoretical calculations show that the fall 
speed can be significantly affected when the aspect ratio 
(ratio of major axis to the minor axis of an ellipsoid) exceeds 
1.5. For example, an aspect ratio change from 1 to 5 reduces 
fall speeds by about 30% and 50% for particles smaller than 
1.0 pm and within the 10-20 pm range, respectively [e.g., 
Cheng et al., 1988; Ginoux, 2003]. Therefore, one would 
expect populations of irregular particles and spherical parti- 
cles of the same mass and density to gradually separate, 
and dust layers to stratify, with the eventual preponderance 
of the most aspherical particles at the top as illustrated 
schematically in Figure 1 . Our observations suggest that such 
processes may actually occur under ambient conditions. 


2. Methodology and Data 

[ 4 ] Here we exploit the fact that the depolarization of lidar 
backscatter signal (echo) is determined primarily by particle 
shape [e.g., Ansmann, 2003]. The depolarization ratio is 
defined as the ratio of an echo’s cross-polarized and 
copolarized components. CALIOP (Cloud-Aerosol Lidar 
with Orthogonal Polarization) provides two types of 
depolarization ratio: volume depolarization ratio and 
particulate depolarization ratio. Since the particulate 
depolarization ratio is free of contributions from molecular 
scattering and attenuation, whereas the volume depolariza- 
tion ratio is influenced by both, the particulate depolariza- 
tion ratio is more suitable for describing the properties of 
dust itself. Therefore, depolarization ratio henceforth refers 
to “particulate depolarization ratio”, denoted as 5. This ratio 
is zero for spherical particles but not for aspherical ones. 
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Figure 1. Schematic illustration of shape-induced diffe- 
rential sedimentation: because of the higher settling speed 
of spherical particles, the initially mixed groups of 
spherical (blue) and aspherical (orange) particles (all at 
the same mass and volume), gradually separate, creating 
a stratified structure with the preponderance of aspherical 
particles at the top. 

[ 5 ] In this study, we use § observations for Saharan dust to 
infer the effects of shape on dust settling and on the resulting 
vertical structures. Specifically, we use dust data provided by 
the CALIOP lidar onboard the Cloud Aerosol-Lidar and 
Infrared Pathfinder Satellite Observations (CALIPSO) 
satellite [e.g., Winker et al., 2003; Vaughan et al., 2004]. 
CALIPSO was launched in 2006 and provides range- 
resolved cloud and aerosol data along its track, including 
total lidar backscatter at 532 nm and 1 064 nm, and perpendic- 
ularly polarized lidar backscatter at 532 nm. CALIPSO 
operational algorithms (currently in Version 3) use this data 


along with altitude and latitude information for feature 
identification and classification, [e.g., Liu et al., 2009; 
Omar et al., 2009]. 

[6] This study is based on nighttime data from the three 
summer months in 2007 over the North Atlantic region 
(0°N-30°N, 0°W-90°W) shown in Figure 2a. This time 
period and region were selected in order to include the typical 
transport path of African desert dust outbreaks. The main 
advantage of nighttime data is a much higher signal-to-noise 
ratio. The study includes about 300 satellite tracks and more 
than 10 6 laser shots within the selected North Atlantic region. 
To detect longitudinal variations in the vertical distribution of 
dust depolarization ratio, we subdivide the North Atlantic 
into nine consecutive strips, labeled sequentially with 
increasing distance from the west coast of Africa 
(Figure 2a). Note that the nine regions are sufficiently large 
to yield adequate statistics with low noise levels, yet 
sufficiently small to constrain regional variability. 

[ 7 ] Let us now examine the vertical distribution of 5 (dust 
particulate depolarization ratio) using CALIOP dust data 
provided at 5 km horizontal resolution. The CALIOP aerosol 
layer product delivers vertically averaged particulate 
depolarization ratio for each aerosol layer by removing the 
contributions from molecular Rayleigh scattering and attenu- 
ation. Our analysis keeps the noise level low by aggregating 
statistics over each of the nine 10° by 30° regions at 270 or 
540 m vertical resolution. A dust sample is counted in an 
altitude bin if the dust layer top and base altitudes in the 
CALIOP aerosol layer product indicate that the altitude 
bin contains dust within a 5 km-resolution segment. The 
depolarization ratio of the dust layer is then included into 
the statistics of this altitude bin. In order to reduce interfer- 
ence from clouds, only dust layers in fully cloud-free 
columns are used (with no clouds identified by either 
the 0.333 km, 1km, or 5 km resolution CALIOP cloud 
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Figure 2. (a) Positioning of the nine consecutive regions studied along the typical summertime track of Saharan dust. Colors 
encode values of the median optical depth of dust layers at 2° by 2° horizontal resolution, (b) Vertical distribution of the 
number of 5 km-resolution dust samples in the nine regions. Each bin has a height of 270 m. 
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Figure 3. Median of dust particulate depolarization ratio 5: (a) as a function of altitude in the nine regions; (b) as a function 
of altitude in Regions 2 to 6; (c) as a function of longitude-regions. The vertical resolution is 540 m. Uncertainties in median 
estimates are calculated by the bootstrapping algorithm [ Efron and Gong, 1983], As the dust leaves the west coast of Africa 
(Region 2), the § values at all altitudes above ~2 km are about the same but in the course of westward transport, the depolar- 
ization spreads out with higher ratios at higher altitudes. Furthermore, the depolarization ratio lapse rate increases westward. 
The altitude-independent § in Region 2 indicates that particles of various shapes are randomly mixed throughout the 
atmospheric column. The progressive westward dispersion of § suggests a separation of dust by shape-induced differential 
settling. The longer time causes wider vertical separation, and consequently, a higher lapse rate of the vertically increasing 
5 is observed with westward transport. See text for more details. 


products within the entire 5 km segment). Finally, we con- 
strain misclassification between ice clouds and dust by consid- 
ering only data below 5 km, as ice clouds during summer in 
the examined region seldom occur below this altitude. 

3. Results and Discussion 

[s] The vertical structure of the number of dust samples 
and of the median particulate depolarization ratio in our 
CALIOP dust data set is shown in Figures 2b and 3. 
Figure 2b indicates that most dust samples are between 
1.5km and 5 km with a median altitude around 3. 0-3. 5 km 
while dust traverses Regions 2 through 5. The dust sample 
vertical distribution changes in Region 6 and most of the dust 
samples are well below 3 km in Regions 7 to 9. In addition, 
more dust samples are found in regions closer to the west 
coast of Africa, except for Region 2 that has a smaller number 
of samples because only half of this region is covered by 
ocean. Because Region 1 is mostly covered by land and the 
number of its dust samples over ocean is very small, this 
region will be excluded from further consideration. It is noted 
that as dust outbreaks occur frequently in summers [ Huang 
et al., 2010], the chosen data set contains many outbreaks 
without having to combine data from different seasons. 

[ 9 ] Figure 3a depicts the vertical profiles of particulate 
depolarization ratio for the dust samples shown in 
Figure 2b. In order to highlight the most important features, 
Figures 3b and 3c present crosssections of the results from 
Figure 3 a for Regions 2-6. The main features can be summa- 
rized as follows. (1) While the number of dust samples peaks 


around 3-3.5 km altitude in Regions 2-6 and below 1 .5 km in 
Regions 7—9, the dust depolarization ratios in all regions 
increase with altitude', (2) In Region 2, closest to west coast 
of Africa, the median depolarization ratio is around ~0.3 
and is fairly constant between 1.5 km and 5 km altitudes. 
This indicates that near the source all shapes are mixed 
throughout the column; (3) In Figures 3b-3c, the dust 
depolarization ratios in each of Regions 2-5 are larger than 
in their eastern neighbors at altitudes above 4 km, but smaller 
at 1 .5-4 km. This implies a clear increase in the lapse rate of 5 
during westward transport. 

[ 10 ] The basic features described above are consistent with 
the shape-induced differential settling of dust particles. To 
illustrate this process, let us consider a minimal model: set- 
tling of only two spheroid particle shapes, with aspect ratios 

and 7 v 2 . The initial vertical distribution of total samples at 
t = 0 time is then f 0 (h)=f° 1 (h)+f°2(h) where f°\(h) and 
f° 2 (h) are the initial vertical distributions of samples for the 
two particle shapes. We further assume that all particles have 
the same density and mass, and hence mass-equivalent size. 
This two-shape, one-size assumption is plausible because 
measurements do not show significant dependence of aspect 
ratio on particle size [e.g., Reid et al., 2003], allowing us to 
disentangle shape and size effects. Because higher aspect 
ratio corresponds to lower settling speed (v) and higher 
depolarization ratio (8), a relationship > ^ 2 in our simple 
dust layer implies vj < v 2 , and 5! > 8 2 . 

[ 11 ] Due to shape-induced differential settling, the two 
kinds of particles, at first evenly mixed throughout atmo- 
spheric columns, will gradually separate during westward 
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Figure 4. Illustration of shape-dependent settling in a minimal two-shape one-size model of differential dust settling, (a) 
Schematic illustration of a shape-induced vertical separation of dust populations, (b) Simulated temporal evolution of depo- 
larization ratio § vertical distribution when the two groups of dust particles are well mixed initially, and have Gaussian vertical 
distributions, (c) The temporal evolution of the depolarization ratio at different altitudes. Figures 4b and 4c display features 
that are similar to those of the CALIOP observations in Figures 3b and 3c. See text for details. Note that this model does not 
include contributions from possible retrieval errors due to instrument/algorithm uncertainties or to mixing or interactions with 
other components in the atmosphere. 


transport, resulting in particle populations with larger aspect 
and depolarization ratios at higher altitudes, as shown in 
Figure 4a. This behavior agrees well with the observations, 
which show larger S values at higher altitudes after dust 
leaves the west coast of Africa (Figure 3). The separation of 
the two groups of particles increases with time, and 
eventually particles of A,i dominate and a maximum 5! is 
reached at high altitudes, while particles of A 2 dominate 
and a minimum 5 2 is reached at low altitudes. Overall, with 
westward progress, 8 increases (decreases) monotonically 
at high ( low) altitudes and the lapse rate of 5 also increases 
as seen in Figure 3. 

[ 12 ] The above scenario of two-shape one-size dust layers 
can be formalized as follows: if the shapes of function/ 0 ^/;) 
and f° 2 {h) do not change significantly during the settling 
process, the vertical distribution of dust samples becomes 
fih, t)=fi(h, t)+f 2 (h, t) where f x {h, t)=f° l {h + vt) and f 2 
( h , t)=f° 2 (h + v 2 »t) are the two dust sample vertical distri- 
butions at time t, correspondingly. Additionally, the aver- 
age 8 can be estimated using 8 {h, t) = [8 i »(l — a) + 82 * 01 ] 
where ct=f 2 (h,t)/[f\(h,t)+f 2 (h,t)\ is the relative weight of 
dust particles with A 2 and ranges from 0 to 1 . 

[ 13 ] Figures 4b and 4c depict the evolution of depolariza- 
tion ratios arising from the shape-dependent settling mecha- 
nism of the minimal model illustrated in Figure 4a. The 
results are qualitatively consistent with those in Figure 3, 
suggesting that the basic features of CALIOP observations 
come from shape-dependent sedimentation. In this numerical 


illustration, we selected parameters based on the following: 
(1) because X most frequently occurs around 1.5 and ranges 
from 1 to 5 [e.g., Ginoux, 2003] and because 8 has relatively 
distinguishable values around A ~1.2 [ Mishchenko and 
Hovenier, 1995], we choose X\=2.1 (at the larger end of 
recent measurements by Kandler et al. [2009] but below 
the average X of 1.2 and 5) and 8j =0.4 [Sakai et al., 2010] 
for the first group of particles, and A 2 = 1 . 1 (the average X 
of 1 and 1.2) and 8 2 = 0.1 for the second group; (2) addition- 
ally, since the typical size is ~7 pm in diameter for transatlan- 
tic dust [Ginoux, 2003], we choose v 2 ~ 1600 cm/hr (nearly 
spherical particle because of A 2 ~l.l) and vi=v 2 (1-0.2) 
where 0.2 is the fall speed reduction rate for Ai and is esti- 
mated from the theoretical reductions in fall speed for A. = 1 
and A = 5 for 7 pm size [e.g., Ginoux, 2003]; and (3) finally, 
the initial distribution functions of f°i(h) and f° 2 (h ) are set 
to be Gaussian functions with a mean altitude of 4.5 km 
and standard deviation of 3 km. The illustrated 8 vertical 
distributions are the averages over 24h-long time spans. 
While a skeptical reader might object that these values 
chosen for the model parameters are somewhat arbitrary, it 
should be noted that the key features of the results hold for 
a wide range of values used in our simulations. The differ- 
ences between the model and CALIOP observations are 
mostly in the specific values such as the altitude of depolari- 
zation ratio at which 8 in a certain region surpasses 8 in the 
preceding region. Generalizing to several aspect ratios and 
size distributions readily reduces these differences. Note that 
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although altitudes below 2 km are displayed in our figures, 
one should not compare observations and simulations below 
2 km, because CALIOP observations are affected by nearly 
spherical marine aerosols at these altitudes. 

[ 14 ] While shape-dependent sedimentation has been 
examined in theoretical studies [e.g., Ginoux, 2003], to our 
knowledge, no link of this process to observations has ever 
been attempted. This may be due to limitations in the 
available observations, such as the large noise in measure- 
ments of 8 or insufficient spatial coverage and duration of 
observations. Our CALIOP data set of transatlantic dust over 
a summer season reveals that the shape-dependent settling 
and the resulting shape-dependent layer structure are signifi- 
cant at oceanic and seasonal scales. This is consistent with 
the earlier studies that characterized African dust during the 
Puerto Rico Dust Experiment and found dust to have a 
significantly higher aspect ratio near 1.9 after long-range 
transport [Reid et al., 2003, p. 7-9], 

[ 15 ] Are there alternative explanations? Consider, for 
example, the formation of loose aggregates that can lead to 
the higher aspect ratios and trigger shape-dependent settling. 
While aggregation could explain the westward increase in 5 
(and may contribute to it), it does not readily explain the 
increase of 5 with altitude. We considered several other 
plausible explanations such as shape softening via humidifi- 
cation, etc. but to no avail. For example, in humidification, 
it seems unlikely that circumstances would conspire to 
supply just the right evolving humidity profiles, with 
westward increase in humidity at higher (but not lower) 
altitudes. Finally, mixing spherical marine aerosols into the 
Saharan Air Layer (SAL) from below is unlikely to cause 
the 8-drop at altitudes between 2 and 4 km in westward 
Regions 2-5 (in Figure 3c), because this mixing at such 
altitudes and regions is rare due to the confinement of marine 
boundary layer by the SAL. 

[ 16 ] Could these CALIPSO observations be spurious, i.e., 
attributed to algorithmic artifacts? Indeed, CALIPSO dust 8 
values are retrieved quantities that mostly rely on a 
predefined lidar ratio. Because this predefined lidar-ratio is 
a constant, it cannot follow the fluctuations of dust size and 
density that play a role in multiple scattering and thus, in 
the ideal lidar-ratio values [Wandinger et al., 2010], As such, 
a constant lidar-ratio could conceivably lead to a westward 
and upward increasing bias of 8s as large particles settle out 
gradually during transport multiple scattering is reduced. 
Then, a higher lidar-ratio is required to compensate for loos- 
ing large particles. However, this mechanism does not lead to 
the westward decreasing trend of depolarization ratio at 2- 
4 km in Regions 2-5 observed by CALIPSO (Figure 3c). 
Moreover, the CALIPSO 8 retrieval is not sensitive to the 
lidar ratio [Liu et al., 2013], and the backscatter coefficient 
retrievals compare well with the ground measurements 
[Wandinger et al., 2010], Therefore, the observed trend in 
particulate depolarization ratio is more likely ascribed to 
shape-induced sedimentation. 

4. Concluding Remarks 

[ 17 ] In this study we asked: Does the effect of particle 
shape on differential sedimentation cause an evolving 
vertical stratification as dust plumes are transported 
westward across the Atlantic? If so, such stratification should 
be discernible from CALIOP data on 8. We then searched for 


and discovered evidences of dust particle asphericity 
affecting the evolution and patterns of scattered light depolar- 
ization ratio (8). 

[is] Specifically, shape-induced gravitational sorting is 
evidenced by the summertime observations of evolving 
patterns in lidar signal depolarization measurements from 
the CALIPSO satellite. Here we demonstrate that even a 
minimal two-shape one-size model of dust particle settling 
captures qualitatively the basic features of observations. 
The results show a clearly increasing stratification in 8 as dust 
moves. Finally, we note that knowledge of dust vertical 
distributions can improve the accuracy of aerosol optical 
depth and particle size retrievals. For example, in optically 
thick dust layers, particle properties near the top can affect 
satellite observations more than those at lower altitudes, thus 
the assumption of vertical homogeneity in dust property 
retrievals can lead to biased results when the vertical distribu- 
tion of particle shape is inhomogeneous. Therefore, our 
results imply that particle shape-dependent settling (hence 
the induced vertical distribution of 8) is important not only 
in studies of dust transport and sedimentation but also in 
estimating dust optical properties, radiative forcing, and 
dust loading. 
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